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We consider the dynamics of a viscous cosmological fluid in the generalized Randall- Sundrum 
model for an anisotropic, Bianchi type I brane. To describe the dissipative effects we use the Israel- 
Hiscock-Stewart full causal thermodynamic theory. By assuming that the matter on the brane obeys 
a linear barotropic equation of state, and the bulk viscous pressure has a power-law dependence on 
the energy density, the general solution of the field equations can be obtained in an exact parametric 
form. The obtained solutions describe generally a non-inflationary brane world. In the large time 
limit the brane Universe isotropizes, ending in an isotropic and homogeneous state. The evolution 
of the temperature and of the comoving entropy of the Universe is also considered, and it is shown 
that due to the viscous dissipative processes a large amount of entropy is created in the early stages 
of evolution of the brane world. 

o 
o . 

(NJ ; I. INTRODUCTION 

o 

D . The idea |LE§ that our four-dimensional Universe might be a three-brane embedded in a higher dimensional space- 
■ time has attracted much attention. According to the brane-world scenario, the physical fields in our four-dimensional 
space-time, which are assumed to arise as fluctuations of branes in string theories, are confined to the three brane, while 
gravity can freely propagate in the bulk space-time, with the gravitational self-couplings not significantly modified. 
The model originated from the study of a single 3-brane embedded in five dimensions, with the bD metric given by 
ds 2 = e^^^ri^ydx^dx" + dy 2 , which can produce a large hierarchy between the scale of particle physics and gravity 
^> due to the appearance of the warp factor. Even if the fifth dimension is uncompactified, standard AD gravity is 
reproduced on the brane. Hence this model allows the presence of large or even infinite non-compact extra dimensions 
^T*" ■ and our brane is identified to a domain wall in a 5-dimensional anti-de Sitter space-time. 

The Randall-Sundrum model was inspired by superstring theory. The ten-dimensional Eg x E% heterotic string 
theory, which contains the standard model of elementary particle, could be a promising candidate for the description 
^vq ' of the real Universe. This theory is connected with an eleven-dimensional theory, M-theory, compactified on the 
orbifold R 10 x S 1 /^ |J. In this model we have two separated ten-dimensional manifolds. For a review of dynamics 
' and geometry of brane Universes see |l2| . 

The static Randall-Sundrum solution has been extended to time-dependent solutions and their cosmological prop- 
erties have been extensively studied ||J|,|,0,|Ul^jn|g|l|,|3 • In 

one of the first cosmological applications of this 
j^Tjl scenario, it was pointed out that a model with a non-compact fifth dimension is potentially viable, while the scenario 
which might solve the hierarchy problem predicts a contracting Universe, leading to a variety of cosmological problems 
. £h ! 101 ' ^ adding cosmological constants to the brane and bulk, the problem of the correct behavior of the Hubble 
parameter on the brane has been solved by Cline, Grojean and Servant jL4|. As a result one also obtains normal 
expansion during nucleosynthesis, but faster than normal expansion in the very early Universe. The creation of a 
spherically symmetric brane-world in AdS bulk has been considered, from a quantum cosmological point of view, with 
the use of the Wheeler-de Witt equation, by Anchordoqui, Nunez and Olsen p5| . 

The effective gravitational field equations on the brane world, in which all the matter forces except gravity are 
confined on the 3-brane in a 5-dimensional space-time with ^-symmetry have been obtained, by using a geometric 
approach, by Shiromizu, Maeda and Sasaki [ I6|JT^ ]. The correct signature for gravity is provided by the brane with 
positive tension. If the bulk space-time is exactly anti-de Sitter, generically the matter on the brane is required to 
be spatially homogeneous. The electric part of the 5-dimensional Weyl tensor Ejj gives the leading order corrections 
to the conventional Einstein equations on the brane. The four-dimensional field equations for the induced metric and 
scalar field on the world- volume of a 3-brane in the five-dimensional bulk with Einstein gravity plus a self-interacting 
scalar field have been derived by Maeda and Wands fjj| . 

The linearized perturbation equations in the generalized Randall-Sundrum model have been obtained, by using the 
covariant nonlinear dynamical equations for the gravitational and matter fields on the brane, by Maartens |19|. A 
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systematic analysis, using dynamical systems techniques, of the qualitative behavior of the Friedmann-Robertson- 
Walker (FRW), Bianchi type I and V cosmological models in the Randall-Sundrum brane world scenario, with matter 
on the brane obeying a barotropic equation of state p = (7 — l)p, has been realized by Campos and Sopuerta p2| , p3| . 
In particular, they constructed the state spaces for these models and discussed what new critical points appear, the 
occurrence of bifurcations and the dynamics of the anisotropy for both a vanishing and non- vanishing Weyl tensor in 
the bulk. 

The general exact solution of the field equations for an anisotropic brane with Bianchi type I and V geometry, 
with perfect fluid and scalar fields as matter sources has been found in [p4f . Expanding Bianchi type I and V brane- 
worlds always isotropize, although there could be intermediate stages in which the anisotropy grows. In spatially 
homogeneous brane world cosmological models the initial singularity is isotropic and hence the initial conditions 
problem is solved [p5| . Consequently, these models do not exhibit Mixmaster or chaotic- like behavior close to the 
initial singularity [ |26| |. Other properties of brane world cosmologies have been considered in |27) and [^8| . 

Realistic brane-world cosmological models require the consideration of more general matter sources to describe the 
evolution and dynamics of the very early Universe. Limits on the initial anisotropy induced by the 5-dimensional 
Kaluza-Klein graviton stresses by using the CMB anisotropies have been obtained by Barrow and Maartens J2^] and 
Leong, Challinor, Maartens and Lasenby |30|. Anisotropic Bianchi type I brane-worlds with a pure magnetic field 
and a perfect fluid have also been analyzed|31fl • Rotational perturbations of brane world cosmological models have 
been studied in fjgfl . 

Most of the investigations of brane cosmological models have been performed under the simplifying assumption of 
a perfect cosmological fluid. But in many cosmological situations an idealized fluid model of matter is inappropriate, 
especially in the case of matter at very high densities and pressures. Such possible situations are the relativistic trans- 
port of photons, mixtures of cosmic elementary particles, evolution of cosmic strings due to their interaction with each 
other and surrounding matter, classical description of the (quantum) particle production phase, interaction between 
matter and radiation, quark and gluon plasma viscosity etc. From a physical point of view the inclusion of dissipative 
terms in the energy-momentum tensor of the cosmological fluid seems to be the best motivated generalization of the 
matter term of the gravitational field equations. 

The first attempts at creating a theory of relativistic dissipative fluids were those of Eckart (3j| and Landau and 
Lifshitz HU . These theories are now known to be pathological in several respects. Regardless of the choice of equation 
of state, all equilibrium states in these theories are unstable and in addition signals may be propagated through the 
fluid at velocities exceeding the speed of light. These problems arise due to the first order nature of the theory, that 
is, it considers only first-order deviations from the equilibrium leading to parabolic differential equations, hence to 
infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle of causality. Conventional 
theory is thus applicable only to phenomena which are quasi-stationary, i.e. slowly varying on space and time scales 
characterized by mean free path and mean collision time. 

A relativistic second-order theory was found by Israel |35| and developed by Israel and Stewart |36| into what is 
called "transient" or "extended" irreversible thermodynamics. In this model deviations from equilibrium (bulk stress, 
heat flow and shear stress) are treated as independent dynamical variables, leading to a total of 14 dynamical fluid 
variables to be determined. However, Hiscock and Lindblom |37) and Hiscock and Salmonson (3£| have shown that 
most versions of the causal second order theories omit certain divergence terms. The truncated causal thermodynamics 
of bulk viscosity leads to pathological behavior in the late Universe, while the solutions of the full causal theory are 
p9| : a) for stable fluid configurations the dissipative signals propagate causally, b) unlike in Eckart-type's theories, 
there is no generic short wave-length secular instability and c) even for rotating fluids, the perturbations have a 
well-posed initial value problem. For general reviews on causal thermodynamics and its role in relativity see f40[| and 

Causal bulk viscous thermodynamics has been extensively used for describing the dynamics and evolution of the 
early Universe or in an astrophysical context. But due to the complicated character of the evolution equations, very 
few exact cosmological solutions of the gravitational field equations are known in the framework of the full causal 
theory. For a homogeneous Universe filled with a full causal viscous fluid source obeying the relation £ ~ p 1 ^ 2 , exact 
general solutions of the field equations have been obtained in [j42j , jl3| , [Q , . In this case the evolution of the bulk 
viscous cosmological model can be reduced to a Painleve-Ince type differential equation, whose invariant form can be 
reduced, by means of non-local transformations, to a linear inhomogeneous ordinary second-order differential equation 
with constant coefficients fiq| . It has also been proposed that causal bulk viscous thermodynamics can model on a 
phenomenological level matter creation in the early Universe |43|] , ||7] . Exact causal viscous cosmologies with £ ~ p s 
have been obtained in [Q . 

Because of technical reasons, most investigations of dissipative causal cosmologies have assumed FRW symmetry 
(i.e. homogeneity and isotropy) or small perturbations around it [49] . The Einstein field equations for homogeneous 
models with dissipative fluids can be decoupled and therefore are reduced to an autonomous system of first order 
ordinary differential equations, which can be analyzed qualitatively ]5Q], pj|. 



2 



The influence of the bulk viscosity of the matter on the brane has been analyzed, for an isotropic flat FRW geometry, 
in [j52| . Dissipative viscous effects in the brane world lead to important differences in the cosmological dynamics, as 
compared to the standard general relativistic cosmology. Since the effects of the extra-dimensions and also the viscous 
effects are more important at high matter densities, the most important contribution to the energy density of the 
matter in this regime comes from the quadratic term in density. Consequently, during the early period of evolution of 
the brane world the energy density of matter is proportional to the Hubble parameter, in opposition to the standard 
general relativistic case with energy density proportional to the square of the Hubble parameter. 

It is the purpose of the present paper to investigate the effects of the bulk viscosity of the cosmological matter fluid 
on the dynamics of an anisotropic, Bianchi type I brane world. In order to solve the field equations, we assume a 
specific equation of state for the bulk viscous pressure. With this choice the general solution of the field equations 
can be expressed in an exact parametric form. The obtained solution describes a non-inflationary Universe, tending 
in a large time limit to an isotropic homogeneous geometry. The behavior of the temperature and comoving entropy 
is also considered. 

The present paper is organized as follows. The field equations on the brane describing the evolution of a viscous 
cosmological fluid are written down in Section II. In Section III we present the general solution of the field equations 
in the case of a power-law dependence of the bulk viscous pressure on the energy density. In Section IV we discuss 
and conclude our results. 



II. DISSIPATIVE COSMOLOGICAL FLUIDS ON THE ANISOTROPIC BRANE 

In the 5D space-time the brane-world is located as Y(X I ) = 0, where X 1 , / = 0,1,2,3,4 are 5-dimensional 
coordinates. The effective action in five dimensions is Jl^] 

t> A l i / j4„, / - / 1 zs± \ i r matter \ 



S = J d^X^ ^i?5 - A 5 J + d^x^—g A± - A + L mattcr j , (1) 

with k 2 = 8irGs the 5-dimensional gravitational coupling constant and where x^ 1 , /i = 0,1,2,3 are the induced 4- 
dimensional brane world coordinates. R5 is the 5D intrinsic curvature in the bulk and is the intrinsic curvature 
on either side of the brane. 

On the 5-dimcnsional space-time (the bulk), with the negative vacuum energy A5 and brane energy-momentum as 
source of the gravitational field, the Einstein field equations are given by 

Gu = k 2 T u , T T j = -A 59IJ + S(Y) [-\ giJ + T?f Uer ] . (2) 

In this space-time a brane is a fixed point of the Z% symmetry. In the following capital Latin indices run in the 
range 0, 4 while Greek indices take the values 0, 3. 

Assuming a metric of the form ds 2 — (njrij + gj^dx 1 dx J , with njdx 1 = d\ the unit normal to the \ = const, 
hypersurfaces and gu the induced metric on \ = const, hypersurfaces, the effective four-dimensional gravitational 
equations on the brane (which are the consequence of the Gauss-Codazzi equations) take the form |l(i| . |r7[ : 

Gftu = Ag^ u + k^T^ v + k^Sftv E^ Vl (3) 

where S^ v is the local quadratic energy-momentum correction 

S„u = ^TT,„ - -T«T va + L 9ta , (3T^T al3 - T 2 ) , (4) 

and E^v is the non-local effect from the bulk free gravitational field, transmitted via the projection of the bulk Weyl 
tensor Ciajb'- Ejj = CiAjBn A n B and Ejj — > E^Sj5j as \ — * 0. The four-dimensional cosmological constant, 

A, and the gravitational coupling constant on the brane, k±, are given by A = fcf (A 5 + ^f-^j /2 and k\ = fc|A/6, 

respectively. 

The Einstein equation in the bulk, the Codazzi equation, also implies the conservation of the energy momentum 
tensor of the matter on the brane, 

A,T/ = 0. (5) 
Moreover, the contracted Bianchi identities imply that the projected Weyl tensor should obey the constraint 
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(6) 



Finally, the equations (|3|)-(|5|) and (|6|) give the complete set of field equations for the brane gravitational field. 

For any matter fields (scalar field, perfect fluids, kinetic gases, dissipative fluids etc.) the general form of the brane 
energy-momentum tensor can be covariantly given as ffO|| 



Tpv = pu^U v + Peffhpy + TT^u + 2q (jl U v y 



(7) 



The decomposition is irreducible for any chosen 4- velocity u M . Here p and p e ff are the energy density and the effective 
isotropic pressure, and = g^ v + u^Uy projects orthogonal to u M . The energy flux and the anisotropic stress 
"K^v obey the conditions = <7< M > and = 7r <M „>, respectively, where angular brackets denote the projected, 
symmetric and trace-free parts: 



V<^> = h»V v 



w. 



<[iv> 



h(jj, a h„y 



.0 



1 



W af3 . 



(8) 



We assume that the equilibrium thermodynamics pressure p of the cosmological fluid obeys a linear barotropic 
equation of state p = (7 — 1) p, 7 G [1, 2] and 7 ^constant. 

The effect of the bulk viscosity of the cosmological fluid on the brane can be considered by adding to the usual 
thermodynamic pressure p the bulk viscous pressure II and formally substituting the pressure terms in the energy- 
momentum tensor by p e ff = P + n. 

The anisotropic stress 7r M „ of the matter on the brane satisfies the evolution equation pTj| 



r 2 t 



(9) 



where rj is the shear viscosity coefficient, r 2 = 2t7/3 2 , with /3 2 the thermodynamic coefficient for the tensor dissipative 
contribution to the entropy density and a a p is the shear tensor. In the following we assume that the main contribution 
to entropy generation is via scalar dissipation, that is, we consider that the dissipative contribution from the shear 
viscosity can be neglected, 77 ~ 0. Consequently, the anisotropic stresses of the matter on the brane also vanish, 
7711/ ~ 0. We consider that in Eq. (0) the heat transfer is zero, that is, we take q^ = 0. 

Then the quadratic corrections to the matter energy momentum tensor on the brane are given by fl6[ | 



Spy = ^p 2 u^u u + -^p(p + 2p e f f )h IMU . 



(10) 



The symmetry properties of E^ v imply that in general we can decompose it irreducibly with respect to a chosen 
4- velocity field as 



II ( Ufj,u u + -h^ ) + V^v + 2(2(^1',, 



(11) 



where k — k 5 /k 4 . The non-local anisotropic fields V^v and Q M include five-dimensional gravitational wave modes. To 
determine the brane dynamics one must solve the 5D field equations completely, by choosing appropriate boundary 
conditions in the bulk. For the only known solution of the field equations in the bulk, the Schwarzschild- anti de Sitter 
bulk that contains a moving brane, the contribution of these terms is zero |J. The effects on the CMB radiation 
anisotropy of the non-local stresses under various plausible physical assumptions have been discussed in Q and Q . 
One such possibility is to assume that the time evolution of the non-local stress is proportional to the energy density 



of the anisotropic source p , with p perturbatively small relative to the matter energy density p, p << p. Since, 
on the other hand, the viscous dissipative effects are proportional to the matter density, their effects will generally 
dominate the non-local effects from the bulk. Therefore we shall also assume a vanishing effective non-local energy 
density, «0« Q M . 

The line element of a Bianchi type I space-time, which generalizes the flat Friedmann- Robertson- Walker metric to 
the anisotropic case, is given by 



ds 2 



dt 2 + a\ {t) dx 2 + a\ (t) dy 2 + a\ (t) dz 2 . 



(12) 



We define the following variables: 



4 



V = f]a i ,3H = y^H i ,H i = — , AHi = H t — H,i = 1, 2, 3. (13) 

i=l i=l 

InEqs. (|l|) V is the volume scale factor, TJ^, i = 1,2,3 are the directional Hubble parameters, and H is the mean 
Hubble parameter. From Eqs. ( |l3| ) we also obtain H = V/3V. 

The physical quantities of observational interest in cosmology are the expansion scalar 9, the mean anisotropy 
parameter A and the shear scalar a 2 , which are defined according to 

9 = 3H, 3A = ]T (^-) 2 , 2<r 2 = a lk a tk = ^ H? - 3H 2 = 3AH 2 . (14) 

»=1 ^ ' i=l 

Using the variables ([l3]), the Einstein gravitational field equations, the Bianchi identity and the evolution equation 
for the non-local dark energy take the form: 

3F + ^^ 2 =A+^[(2-3 7 )p-3n] + ^[(l-3 7 )p-3n]-|i, (15) 

2—1 * 

~ t (VH t ) =A+|[(2- 7 )p-n] + ^[(l- 7 )p-n] + |^,z^l,2,3, (16) 



p + 3-yHp = -3HTL, (17) 
U + AHU = Q. (18) 

The causal evolution equation for the bulk viscous pressure II is given by [ [i"H 

T fl + U = -3^H-^rU^3H+^-^-^j, (19) 

where r is the relaxation time, T is the temperature and £ is the bulk viscosity coefficient. 

Eq. (|l9| ) arises as the simplest way (linear in n) to satisfy the i7-theorem (i.e., for the entropy production to be 
non-negative, S^- tli = II 2 /(£T) > 0). The particle flow vector is given by = nu^, where n > is the particle 
number density. In the framework of causal thermodynamics and limiting ourselves to second-order deviations from 
equilibrium, the entropy flow vector takes the form — eN^ — riI 2 w AI /2£T, where e is the entropy per particle. 

An important observational quantity is the deceleration parameter 

d ( 1 \ p + 3 P + 3n + 2^ 
9=37 77 - 1 = Tt . i<2, ■ (20) 



The sign of the deceleration parameter indicates whether the model inflates or not. The positive sign of q corresponds 
to "standard" decelerating models, whereas the negative sign indicates inflation. 



III. GENERAL SOLUTION OF THE FIELD EQUATIONS 

Integrating dlg| ) yields for the dark matter term the expression U = UqV~^ , where Uq > is an arbitrary constant 
of integration. 

By adding Eqs.(|l6|) we find 

7^(^)=^ 2 = A+f [(2-7)p-H] + f [d-7)p-n ] + ^. (21) 



Substracting Eq. ( |2l| ) from Eqs. ( J16| ) we obtain 

H i = H+^,i = 1,2,3, (22) 

with Ki, i = 1, 2, 3 constants of integration satisfying the consistency condition X)»=i Ki = 0- By using Eq. (p2|), it is 
easy to show that 
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I>? = 3tf 2 + fJ, (23) 

i=i 

where we denoted K 2 = J2i=i K f- 

Usually, in order to find the general solution of the gravitational field equations, the equations of state of the bulk 
viscosity coefficient, temperature and relaxation time are given, in terms of the energy density of the cosmological 
fluid. Then the thermodynamic and geometric parameters are obtained by solving the field equations. But in the 
present paper we shall use a different approach, namely we shall assume that the functional dependence of the bulk 
viscous pressure on the energy density is known. Therefore the behavior of the bulk viscosity coefficient, temperature 
and relaxation time can be obtained from the field equations. Hence for the bulk viscous pressure we assume the 
following general power-law type density dependence: 

If = — 7ir 0/ o m , to — constant > 0, m ^ 1, Ho = constant > 0. (24) 

This form of the bulk viscous pressure is sugested by the form of II in the non-causal theory, where II ~ £H . This 



relation follows by taking r = in Eq.(19). Assuming a power law density dependence of the bulk viscosity coefficient 
on the density, £ ~ p s , s =constant> 0, and taking into account that generally the Hubble parameter can also be 
approximated by some power of the energy density, we obtain Eq. ( p4| ) in a natural way. By using such a form for the 
viscous stress Barrow | p3| has obtained a class of homogeneous and isotropic cosmological models which begin in a 
de Sitter state but subsequently deflate towards the flat Friedman model. Barrow Q has also considered the effects 
of the anisotropics on the early cosmological evolution, within the framework of the non-causal thermodynamical 
approch, for two distinct classes of models corresponding to the choices m > 1 and to < 1 in the expression of the 
bulk viscous pressure. In both cases the initial anisotropies are smoothed out by the inflationary expansion of the 
Universe. The range of values of to can be roughly estimated by considering the isotropic limit of the model. Since 
for a high density cosmological fluid s < 1/2 |5J|, and for standard cosmological models H ~ p 1 ^ 2 , it follows that in 
conventional general relativity to < 1. However, in the brane world model, at high densities the quadratic correction 
term dominates the cosmological evolution and the Hubble parameter is proportional to the energy density, H ~ p. 
Therefore values of m greater than one, to > 1, are also allowed to describe high density bulk viscous cosmological 
fluids. For low density viscous matter s — ► 1 J55| and consequently on the brane to could take values as high as 2, 
m -> 2. 

The Israel-Stewart theory is derived under the assumption that the thermodynamic state of the fluid is close to 
equilibrium, which means that the non-equilibrium bulk viscous pressure should be small when compared to the local 
equilibrium pressure, that is |H| <C p. Consequently, the energy density of the cosmological fluid must satisfy for all 
times the condition p « [(7 — 1) /37Ho] 1 ^ m ^ . 



From Eq. fll7|), it follows that the time evolution of the energy density of the matter on the brane is given by 

r - » il/(l-m) 

p(t)=[p v^ m -v + n \ , (25) 

where po > is a constant of integration. 

In the case to = 1 the Bianchi identity takes the form p + 3j e f fH p = 0, where J e ff = 7(1 — Ho). Hence for Bianchi 
type I brane world filled with a viscous dissipative fluid with bulk viscous pressure proportional to the energy density, 
the dynamics is described by the same equations as for the perfect fluid case |2j], with 7 — * j e ff- 



Substituting U =UqV~3 into Eqs. ( |15| ) and ( |21| ) and eliminating n from Eqs. ([15]) and (21), we obtain the 
basic equation for V describing the dynamics and evolution of the Bianchi type I brane world in the framework of 
Israel-Stewart-Hiscock causal thermodynamics: 



K 2 „ A „, a k 2 2 12W 



6H 2 - — = 2A + 2kip+^-p 2 + -^V~3. (26) 



In view of the definition of the mean Hubble function H and Eq. (^5|) , Eq. ( p6| ) can be integrated to give 



1 - i" = \l \ x 



J K 2 + 2AV 2 + 2k 2 4 V 2 (poV*™- 1 ) + Tic) + ^V 2 (poV**™- 1 ) + n f 
with to an arbitrary constant of integration. 



V 3 



dV, (27) 



G 



The scale factors <Zj, i — 1,2,3, are obtained from Eq. ([22|), and are given by 

Oi(t) = a i0 vh Ki I v? dV ,i = 1, 2, 3, 

where a,io, i = 1,2,3 are arbitrary constants of integration. 

The anisotropy parameter can be expressed in terms of V and H in the form 



(28) 



K 2 + 2AV 2 + 2k 2 4 V 2 (p Q V^ m -V + n ) llm + ^V 2 (p Q V^ m -V 



no) 1 - 



12Z4 , 
fcfA 



The expansion scalar 9 can be represented, in a parametric form, as a function of V, as 

~k 2 + 2Av 2 + 2k 2 v 2 (poV^-v + n ) 1/(1 ~ m) + f v 2 ( Po v^ m -v + n ) 2/(1 ~ m) + 



1/2 



The shear scalar a 2 is given by 



V 



K 2 
2V 2 ' 



As a function of V the deceleration parameter q can be expressed as 

(37-2) (poV^-V +Tl ) mi - m) 

q = 



37110//* 



2^ 



(p l/7(m-l) + n ) 



l/(l-m) jf2_ 
2V 2 



(29) 



(30) 



(31) 



(32) 



Integrating Eq. (19), we obtain the time dependence of the temperature of the matter on the brane in the form 



T(t) = T ^!exp 



1 P 

t 7-n 



3iJ . , 

dt 

2 



(33) 



with To an arbitrary constant of integration. 

Following Belinskii, Nikomarov and Khalatnikov |5^| , we suppose that, in order to guarantee that the propagation 
velocity of bulk viscous perturbations, i.e., the non-adiabatic contribution to the speed of sound q, does not exceed 
the speed of light, the relation r = £/c 2 (p + P) — £./l c bP holds in the dissipative fluid. With this choice for the 
relaxation time r the causal evolution Eq. (B3l) gives the following evolution law for the temperature: 



T{t) = T o7 ^p- 



TT 2 

LL „2m-l 



exp 



IP 

7Cfc7 



3^V- 



H 



3H 



dt 



(34) 



An analysis of the relativistic kinetic equation for some simple cases given by Murphy [|56[ , Belinskii and Khalatnikov 
p7| and Belinskii, Nikomarov and Khalatnikov |35| has shown that in the asymptotic regions of small and large values 
of the energy density, the viscosity coefficients can be approximated by power functions of the energy density with 
definite requirements on the exponents of these functions. For small values of the energy density it is reasonable to 
consider large exponents, equal in the extreme case to one. For large p the power of the bulk viscosity coefficient 
should be considered smaller (or equal) to 1/2. Hence, we shall assume that the bulk viscosity coefficient obeys the 
simple phcnomenological law 



(35) 



with a > and < s < 1/2 constants. 

Therefore we obtain for the temperature of the dissipative viscous cosmological fluid on the brane the following 
representation: 



T = T 7^p 2m - 1 exp 

4 



a 



n 2 

il „2m-l 



exp 



Cv ' In' z 



-dV 



K' 2 + 2AV 2 + 2k 2 V 2 (p Vr(™-V + n ) ^ + § V 2 (poV^™- 1 ) + U ) lSm + 



(36) 



7 



The growth of the total comoving entropy £(t) over a proper time interval (to,t) is given by 
where fcs is the Boltzmann constant. 



(37) 



IV. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the evolution of a causal viscous dissipative cosmological fluid in the brane 
world scenario. As only source of dissipation we have considered the bulk viscosity of the matter on the brane. In this 
case the standard Einstein equations are modified, due to the presence of the terms from extra dimensions, quadratic 
in the energy density. 

By assuming a power-law dependence of the bulk viscous pressure on the energy density of the cosmological fluid, 
the general solution of the field equations for an anisotropic Bianchi type I geometry can be obtained in an exact 
parametric form. By assuming for the bulk viscosity coefficient and for the relaxation time the usual equations of 
state, the temperature evolution of the cosmological fluid can also be obtained in a closed form. 

The cosmological evolution of the bulk viscous brane Universe is expansionary, with all the scale factors monoton- 
ically increasing functions of time. For t — > oo we have a» = a = aoV 1 / 3 , i = 1,2, 3. Hence the Universe ends in an 
isotropic state. The expansion parameter 9 is a decreasing function of time, with a singular behavior at t = 0. 

The time evolution of the energy density of the cosmological fluid essentially depends on the value of the exponent 
m in the bulk viscous pressure-energy density relation. The variation with respect the cosmological time of p is 
represented, for m > 1, in Fig. 1. 
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FIG. 1. Time evolution of the energy density p(t) of the Zeldovich fluid (7 = 2) filled Bianchi type I brane Universe, for 
different values of the parameter m: m = 1.25 (solid curve), m = 1.45 (dotted curve) and m = 1.6 (dashed curve). We have 
normalized the constants so that K 2 = 2, A = 0.5, po = Ho, 2fc|rfto ( 1 /' 1 m) = 1, k\p\l^ m ' = A and 12Wo = k\\. 

For this range of values of m the Universe starts its evolution with finite values of the energy density, p(0) = 

Hq^ 1 m ' 7^ 0. But the scale factors are singular at the initial time, ai(t a ) = 0, i = 1, 2, 3. For m < 1 the anisotropic 
brane Universe starts from a singular state, with infinite energy density and zero scale factors. For all values of m 
the energy density is a monotonically decreasing function of time. 

The evolution of the deceleration parameter q is represented in Fig. 2. 
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FIG. 2. Dynamics of the deceleration parameter q(t) of the Zeldovich fluid (7 = 2) filled Bianchi type I brane Universe, for 
different values of the parameter m: m — 1.3 (solid curve), m = 1.4 (dotted curve) and m = 1.5 (dashed curve). We have 
normalized the constant parameters so that K 2 — 2, A — 0.5, po = IIo, 2k%pg ~ m = 1, k\p^ = A and 12Wo = k 2 \. 

Due to the smallness of the bulk viscous pressure, satisfying the condition |IT| « p, the cosmological evolution is 
non-inflationary, with q > for all values of m and for all t. In order to have inflationary evolution it is necessary 
that the negative bulk viscous pressure dominates the thermodynamic pressure. 

The dynamics of the mean anisotropy parameter A strongly depends on the value of m. The time variation of A 
for m < 1 is represented in Fig. 3. In this case the Bianchi type I Universe starts its evolution from an isotropic (but 
singular) state, with A(0) = 0. For small times the mean anisotropy is increasing, reaching a maximum at a finite 
time t = t c . For times t > t c the anisotropy is a monotonically decreasing function of time, and thus in the large 
time limit the anisotropic brane world Universe will end in an homogeneous and isotropic state. As has already been 
pointed out in |24f| , this behavior is specific for brane world cosmological models and cannot be found in conventional 
general relativistic scenarios. For m > 1, the brane Universe starts from a state of maximum anisotropy, also reaching 
in the long time limit an isotropic phase. Due to the expansion of the Universe, the shear scalar a 1 vanishes in the 
large time limit, a 1 — > for V — > 00. 
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FIG. 3. Evolution of the mean anisotropy parameter A(t) of the Zeldovich fluid (7 = 2) filled Bianchi type I brane Universe, 
for different values of the parameter m: m = 0.7 (solid curve), m = 0.8 (dotted curve) and m = 0.9 (dashed curve). We have 
normalized the constants so that K 2 = 2, A = 0.5, po = IIo, 2k 2 , p^ 1 m ' = 1, klp 2 /^ 1 m ' = A and 12Wo = k 2 X. 

According to the formal definition of the isotropization given by Collins and Hawking j58j a cosmological model 
approaches isotropy if the following four conditions hold as t — > 00: i) the Universe is expanding indefinetely and 
H > ii) T 00 > and T 0l /T 0Q — * 0. T oi /T QO represents an average velocity of the matter relative to the surfaces of 
homogeneity. If this does not tend to zero, the Universe would not appear homogeneous or isotropic iii) the anisotropy 
in the locally measured Hubble constant a/H tends to zero, a/H — > and iv) the distortion part of the metric tends 
to a constant. All these conditions are satisfied in the present model, since in the large time limit the bulk viscous 









brane Universe is expandin g an d the energy density of the matter is positive for all times. Moreover, as one can 
immediately see from Eq. ( |l4|) the quantity a/H is proportional to the square root of the anisotropy parameter, 
a/H ~ \f~A and in the large time for all values of the parameter m it tends to zero, a/H — ► 0. As for the condition 
iv), it is also satisfied since in the metric we have considered only the volume part. Spatially homogeneous models 
can be divided in three classes: those which have less than the escape velocity (i.e., those whose rate of expansion is 
insufficient to prevent them from recollapsing) , those which have just the escape velocity and those which have more 
than the escape velocity J58| . Models of the third class do not tend, generally, to isotropy. In fact the only types 
which can tend toward isotropy at arbitrarily large times are types /, V, VIIq and Vllh- For type Vllh there is no 
nonzero measure set of these models which tends to isotropy [p8| . In this sense the Bianchi type I model we have 
studied in the present paper is a very special case. The Bianchi types that drive flat and open Universes away from 
isotropy in the Collins-Hawking sense are those of type VII. 

As a consequence of its energy density dependence, the bulk viscous pressure is a decreasing function of time for 
all values of the parameters and for large times it tends to zero. The time evolution of n essentially depends on the 
constant m. For m >> 1, the time scale in which n « is very short. 

The variation of the temperature, presented in Fig. 4, shows that at the initial times of the cosmological evolution 
the temperature is an increasing function of time, reaching a maximum value at a finite time t c . For times t > t c , the 
temperature is a decreasing function of time. 
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FIG. 4. Time evolution of the temperature T(t) of the Zeldovich fluid (7 = 2) filled Bianchi type I brane Universe, for 
s — 1/4 and for different values of the parameter m: m = 0.45 (solid curve), m = 0.65 (dotted curve) and m = 0.9 (dashed 
curve). We have normalized the constants so that K 2 = 2, A = 0.5, po = IIo, 2k 2 pQ^ 1 m ' = 1, klp 2 /^ 1 = A, VZUq = k 2 X, 
7T0IT0 = <? b , 7Cft = a and 3c 2 = U . 

The comoving entropy, represented in Fig. 5, is an increasing function of time in the early stages of cosmological 
evolution, but it tends to a constant value in the large time limit. 
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FIG. 5. Variation, as a function of time, of the comoving entropy E(t) of the Zeldovich fluid (7 = 2) filled Bianchi type I 
brane Universe, for s — 1/4 and for different values of the parameter m: m — 0.45 (solid curve), m = 0.65 (dotted curve) and 
m = 0.9 (dashed curve). We have normalized the constants so that K 2 = 2, A = 0.5, po = IIo, 2k\p\j^ m ' = 1, klp 2 /^ 1 m ' = A, 
12Wo = fcf A, 7T0II0 = cl, 7Cj = a and 3c^ = IIo • 



The transition from the anisotropic to the isotropic state of the brane world is associated to a large increase in the 
entropy. The increase of the comoving entropy of the cosmological fluid is independent of the value of the parameter 
m and is a general feature of the model. In the large time limit II — > and, as one can see from Eq. (|37]) , in this 
limit £ — > const. 

The equation of state of the Zeldovich fluid with 7 = 2 is equivalent to the equation of state of a massless scalar 
field (p, with p = p = <fi 2 /2 and a dissipative coupling J53| , with 4> satisfying the equation 

4> + 3H^3(2 1 - m )Uo4> 2m - 1 H. (38) 

The evolution of the scalar field can be obtained again in parametric form. Classically, bulk viscous stresses arise 
because the expansion of the Universe is continually trying to pull the fluid out of thermal equilibrium, and the 
cosmological fluid is trying to relax back. Hence the viscosity of a classical fluid arises from the differential motion of 
fluid elements. But in the collisionless periods after the Planck era and during inflation there are no particle interactions 
responsible for the viscosity. However, the phenomenological description of quantum particle production processes 



can be consistently modeled by means of bulk viscosities [|47|, |59|) . In particular, the full causal thermodynamics can 
be reformulated as a theory describing particle production [ 47|. In this formulation the causal bulk viscous pressure 
n in the energy conservation equation ( [j"7| ) acts as a creation pressure p c , p c — H, while the particle production rate T 
is proportional to the bulk viscous pressure divided by the locally measured Hubble parameter. The particle balance 
equation, describing the variation of the particle number due to the combined effects of the expansion of the Universe 
and particle production, takes the form 

n + 3Hn = Tn, (39) 

where T = -U/-/H @. 

As an example of particle creation processes on the brane we consider particle production from a scalar field obeying 
an equation of state of the form p^ = (7^ — 1) p^, < 74, < 1, where p^ = <f> 2 /2 + U (4>) and p^ = (j) 2 /2 — U (</>), with 
U (</>) the self-interaction potential. The deviation from the conformal invariance of the field or the polarization (trace 
anomaly) of quantized fields interacting with dynamic spacetimes will lead to particle production | ^9[ . The rate of 
particle creation in a particular mode depends on its frequency in relation to the expansion rate of the Universe: at 
any given time t in the quantum regime, new particles with frequency lu < t^ 1 are created from the quantum vacuum, 
while particles produced earlier interact and have an expansionary dynamics. Whether the newly created particles 
can reach effective thermal equilibrium depends on the rate of their interaction relative to the rate of expansion. 

To obtain a phenomenological classical description of the essentially quantum creation process, we assume again 
that the negative creation pressure is proportional to the energy density of the scalar field, p c = —PcQl^P™ ', where, 
for simplicity, we suppose that all p c g, 70 and m are constants. Then the time variation of the scalar field is given by 

r 1 1/(1— m) 

P<t>{t) = [Poef'V* +p c0 J , (40) 

where pQ<p > is a constant of integration. The rate at which particles are created on the brane is 

T(t)= Pca H- 1 po<f,V^ m -V+ Pc0 , (41) 



while the particle number varies as 



n(t) = ^exp (J T(t)dt) 



(42) 



with no > a constant of integration. At high temperature n = g( (3) T 3 /ir 2 , where £ (s) is the Riemann function, 
C ( s ) = SfeLi k~ s , with C (3) = 1.202 and g is the total spin state of the particles. Therefore the temperature of the 
very early Universe is given by 

T(t) 1/3 ggfl/EWg) m) 
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In the particular case 7^ = 0, corresponding to a dissipative scalar field obeying an equation of state of the form 
P4> +P4> — l'O, the time variation of the energy density of the field is ~ t 1 /( 1_m ) j while the creation pressure varies 
as p c ~ ^m/f 1 -" 1 ) , j n order to have a decaying field it is necessary that m > 1. Even in the presence of particle creation 
processes, the initial expansion of the brane Universe is non-inflationary. But in the large time limit, for large values 
of V and for a decaying creation pressure with p c — > 0, the deceleration parameter becomes q — [(37^ — 2)] /2. For 
70 < 2/3, the brane Universe ends in an inflationary epoch. 

Since the dynamics on the brane is different from the standard general relativistic one, the particle creation processes 
are also influenced by the effects of the extra dimensions. Because the rate of expansion is higher, the particle creation 
processes are also accelerated in the models with large extra dimensions. 

In the present paper we have pointed out an other important difference between brane world cosmology and standard 
general relativity predictions about the very early Universe, namely the behavior of the temperature of a realistic 
cosmological fluid. In brane world cosmological models, the high density Universe starts from a low temperature state. 
In the initial stages of expansionary evolution the increase of the degree of anisotropy is associated to an increase of 
the temperature of the Universe, the cosmological fluid experiencing a heating process. An other interesting feature 
of the model is the presence of a non-singular energy density, associated to a singularity of the scale factors. The 
present approach can also lead to a better understanding of the entropy generation mechanisms in the very early 
Universe, where viscous type physical processes have probably played an important role even when the Universe was 
about 1000s old II . 

Hence dissipative viscous brane world cosmological model with an anisotropic geometry could maybe describe a 
well-determined period of the very early evolution of our Universe. 
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